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We derive the frequency-resolved heat current expression in the linear response regime for a setup 
comprised of reservoir, interacting central site, and tunneling barrier under the action of a time de¬ 
pendent electrical signal. We exploit the frequency parity properties of response functions to obtain 
the heat current expression for interacting quantum conductors. Importantly, the corresponding 
heat formula, valid for arbitrary AC frequencies, can describe photon-assisted heat transport. In 
particular, we analyze the heat transfer for an interacting multilevel conductor (a carbon nanotube 
quantum dot) coupled to a single reservoir. We show that the electrothermal admittance can reverse 
its sign by properly tunning the AC frequency. 


I. INTRODUCTION 


The expeditious advance of circuit miniaturization re¬ 
quires the knowledge of heat flow in quantum devices 
in response to electrical fieldsji Due to this necessity 
thermoelectrical transport in quantum devices is a re¬ 
search field that has vigorously relaunched nowadays.— 
Nevertheless, so far, the activity has focused on driv¬ 
ing electrical currents by oscillatory forces^i^ for charge 
qubit manipulation^, quantum emitter generation^ or 
quantum tomography purposes.— Low frequency mea¬ 
surements of the electrical admittance in quantum RC 
circuits provide information about the spectroscopic {Cg 
quantum capacitance) and resistive {Rg relaxation re¬ 
sistance) properties of coherent conductors'^— in which 
Rg = /i/2e^ takes an universal value— r2£ Interactions, 
such as charging effects^^ or Kondo correlationsi^iiS 
do not alter the universality of Rg, measured ex¬ 
perimentally, in quantum capacitors carbon 

nanotubeS)^ superconducting junctions;^ and quantum 
dots.^ 

In contrast, time-resolved heat transport has been 
poorly investigated— Solely, stationary or time- 
averaged heat flows have been analyzed in more detail.— 
The understanding of time dependent heat currents 
opens an avenue of creating circuit architectures where 
heat absorption or emission events are finely tunable via 
electrical and thermal time dependent signals. Recently, 
the linear response for the charge and heat fluxes to 
electrical and thermal AC signals was computed for a 
quantum capacitor showing that heat flows can be de¬ 
layed or elapsed with respect to the AC pulse depending 
on the dot gate position.— Later on, Ludovico et ali^ 
showed, for a slow AC modulation and noninteracting 
conductors, that time-dependent heat flow jR{t) needs 
to explicitly consider the energy stored and relaxed at 
the tunnel coupling region when a tunnel Hamiltonian 
description is employed. The goal that we face consists 
of the calculation of the heat flux in a completely dif¬ 
ferent regime. We derive the time dependent heat flow 
for interacting and multi-orbital conductors and arbi¬ 



FIG. 1: Schematic of the multi-orbital interacting quantum 
capacitor. It is shown two situations where the AC voltage 
is applied either (a) to the dot or (b) to the reservoir. Both 
cases are equivalent by means of a gauge transformation. U(t) 
denotes the time-dependent internal potential of the dot as 
response of the charge injected by the AC driving field V{t). 
We model such response with a capacitance C. We illustrate 
in (b) the induced photon-assisted emission and absorption 
processes (shown by wiggle lines) in the heat by the action of 
y{t). Qr, Qd, and Qt are the energy change rates in time 
at the reservoir, quantum dot, and tunneling barriers. 


trary AC frequencies. Since our heat formula holds for 
arbitrary high frequencies, it is able to describe photon- 
assisted tunneling processes occurring in the heat trans¬ 
port. Our only limitation is that we are restricted to the 
linear response regime, i.e., low AC driving amplitudes. 
We confine our interests to the Coulomb blockade regime 
and do not consider higher-order correlations like Kondo 
effect, we implicitly assume a temperature higher than 
the Kondo scale (Tr), i.e., T ^ Tr. For our calcula¬ 
tions, we employ the nonequilibrium (Keldysh) Green’s 
function formalism^S which allows us to include electron- 
electron interactions, charging effects, in a feasible way 
within the so-called Hartree-Fock regime. 

To address these issues, we focus on a prototypical 
interacting multi-orbital quantum circuit- a quantum 
capacitor- formed by a carbon nanotube quantum dot 
that is coupled to a single terminal being modulated by 
an AC voltage, as shown in Fig. [TJ Carbon nanotubes 
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exhibit charging effects due to the formation of quantum 
dots inside the tubei^ The valley degrees of freedom, 
corresponding to the K and K’ Dirac points in graphene 
in addition to the spin indices, lead to a four-fold en¬ 
ergy level degeneracy. Such degeneracy can be lifted 
by the presence of an external magnetic field B.— Fur¬ 
thermore, the nanotube curvature yields a spin-orbit in¬ 
teraction resulting in split time-reversal dot level pairs. 
Therefore, carbon nanotube quantum dots act as multi¬ 
orbital interacting conductors where Coulomb interac¬ 
tions give rise to Coulomb blockade phenomena. Even 
more importantly, carbon nanotube quantum dots have 
been demonstrated to be perfect platforms to investigate 
the frequency-resolved transport when they are embed¬ 
ded in electromagnetic cavities,— precisely the issue that 
our work addresses. 

In Sec. m we describe the theoretical model for 
an interacting multi-orbital quantum capacitor and the 
Hartree-Fock decomposition of electron-electron interac¬ 
tion. We derive the time-dependent heat flux for the 
reservoir in the linear response regime and obtain the 
electrothermal admittance. For completeness we also dis¬ 
cuss the electrical admittance. In Sec. nil we start by 
proving that our formalism can exhibit Coulomb block¬ 
ade phenomena. The observations for a single orbital 
quantum dot are briefly discussed and then results for a 
multi-orbital one are presented. Here, we consider a four 
level degenerated carbon nanotube quantum dot which 
splits due to the presence of magnetic field and spin-orbit 
interaction. For this system, we show the RC parameters 
and investigate the electrical and electrothermal admit¬ 
tances. Our findings are summarized in Sec. IIVI 


II. TIME-DEPENDENT TRANSPORT 
FORMULATION 

The purpose is to derive the linear response heat cur¬ 
rent for an interacting conductor coupled to a reservoir 
that oscillates with AC voltage signal in the Coulomb 
blockade regime. We first propose a model Hamiltonian 
describing a multi-orbital quantum capacitor which is at¬ 
tached to a single reservoir. Using the Keldysh formalism 
for nonequilibrium Green’s functions, we have the linear 
response heat flux at the reservoir subject to an oscillat¬ 
ing electrical signal. Remarkably, the heat flux can be 
cast in terms of the conductor Green’s function. 


A. Model of Hamiltonian 

The starting point of our derivation is the tunnel 
Hamiltonian description of a quantum capacitor. The 
quantum capacitor Hamiltonian R is split into three 
parts, namely the reservoir part {Rr), the dot contri¬ 
bution {Rd), and the tunneling term (Rt), be.. 


More concretely, the reservoir part reads 

'^R = Y1 - (m + eH(t))] cl^Cka , (2) 

k,(T 

where ^ = ErR eVdc is the chemical potential, e > 0 the 
unit charge, and V{t) is the electrical voltage modula¬ 
tion [see Fig.[T](b)]. As shown in Fig. lU notice that this 
description is equivalent to the situation where the AC 
voltage is applied to the dot since a gauge transformation 
connects both situations and we thus employ two situa¬ 
tions interchangeably. Without loss of generality, we set 
Vdc = 0. The operator cl^{cka) creates (annihilates) an 
electron with wavevector k and spin a in the reservoir. 
For the dot contribution describing an interacting system 
with n levels, we have 

Rd = ^ inrjd\„dncT + Ec [Nd + A/’(t)]^ , (3) 

n,tT 

where dl^^{dna) corresponds to the creation (annihila¬ 
tion) operator for a dot electron in the nth level with 
spin cr, Cna- denotes the single-particle energies, and 
Ec = I2C is the electrostatic charging energy (C is 
the dot geometrical capacitance). Here, the dot occu¬ 
pation operator reads Nd = Thn (jdliadnrr- Finally, the 
tunneling term hybridizes the reservoir and the dot sub¬ 
systems according to 

Rt = E inker + i?.C^ . (4) 

n,k,<T 

where tnka denotes the tunneling amplitude. 

When charge is injected in the dot by an external 
source V{t), a polarization charge is created to keep the 
dot as a neutral charge object. In a simple electrostatic 
picture, we model such polarization charge with a capac¬ 
itance C leading to^^ 

eMd{t) = CU{t) , (5) 

in which U{t) is the internal potential of the dot giving 
rise to a time-dependent potential inside the dot. Then, 
the dot Hamiltonian in Eq. up to linear order in 

U{t), can be written as 

Td-D ='^ [ina + eU{t)] dhadner + EcNi. ( 6 ) 

n,(7 

Since Aj = En,a' dl.^dmad'l^rdna' is a quartic op¬ 
erator, the Hamiltonian cannot be solved without intro¬ 
ducing some proper approximation. 

B. Hartree Fock approximation 

We thus perform the Hartree-Fock approximation in 
the dot Hamiltonian. For the Hartree approximation, 
AJ is decoupled in the form 

Hartree — 2 EE d]riadma{d}^iyidna') i (7) 

m,cT n,(T' 


R = Rr + Rd + Rt ■ 


( 1 ) 
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while for the Fock approximation one has 

[-^dlpock = ^ ^ d^m^dna' {dl^^dna') ■ (8) 

m,(7 n,(7' 

Considering Hartree and Fock approximations (HF) to¬ 
gether, the dot Hamiltonian then becomes 

Hd = E ^mntT n<T 1 (9) 

where 

^mncr(^) — I emo- + eC/ (t) 2Ec E 

\ j 

— 2Ec (dl^i^dncr), (10) 

with CT =), / t for a =t / i- 

An important observation is in order when Eq. (0 
is considered. For illustration, we assume that there is 
only a single orbital {cma = ^da)- We allow the pres¬ 
ence of a small external magnetic field eda- = Cd + crA^ 
{Az the Zeeman energy) in order to break explicitly spin 
degeneracy, i.e., Nda ^ Ndg. For a single orbital, since 
the exchange interaction is absent between electrons with 
same spins the Fock term disappears such that the dot’s 
energy level can be simplified as 

eda(i) = Erf + crA^ -t eU{t) + 2Ec{d\ds). (11) 

We envision the level is occupied by the spin cr electron 
with energy €da + eU{t). When another electron with 
opposite spin a enters the level, its energy is increased 
by the charging energy 2Ec- Therefore, as a function 
of the Fermi energy Ep the dot occupation then does 
not change continuously, but shows plateaus and discon¬ 
tinuous jumps. These jumps are the clear evidence of 
Coulomb blockade phenomenon. We illustrate how HF 
approximation captures charging effects by considering a 
single orbital quantum dot and plotting its total occu¬ 
pation Nd when the Fermi energy varies in Fig. [2] We 
consider the absence of AC signal. The dot occupation 
shows a plateau of ~ 1 in the Coulomb gap region, 
i.e., €d < Ep ^ ed + 2Ec- This result corroborates the 
fact that our HF description can reproduce charging ef¬ 
fects and thus Coulomb blockade properly. We emphasize 
that our results are not equivalent or expected when they 
are compared with previous findings using noninteracting 
models. 


C. Heat admittance 


Our main interest is focused on the rate of heat at 
the reservoir. For such purpose, we evaluate the time 
derivative of each component of the Hamiltonian which 
is given by 
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FIG. 2: Quantum dot occupation Nd versus the Fermi energy 
Ep for a single orbital dot in the Coulomb blockade regime. 
The dot’s energy level is placed at tda = Cd + nAz with 
Ed = 0 and Az = O.SF. Rest of parameters: 2Ec = lOF 
and fcflT = O.OdF. 


where S = {R,D,T} is referred to R reservoir, D dot, 
and T tunneling term. Notice that in Eq. m the last 
term is the power supplied by an external source and 
must be subtracted from the definition of the heat rate 
Qs{t) (see below). Besides, since the Hamiltonian oper¬ 
ator H commutes with itself it is fulfilled 

^{[n,n]) = QR + QD + QT = o. (13) 

As mentioned, our goal is to compute Qr. When 
Coulomb interaction is taken into account, we obtain the 
rate of heat at the reservoir, Qr, in terms of interact¬ 
ing dot Green’s functions. Since the direct calculation of 
Qr is cumbersome, our strategy is first to compute the 
energy change rates at the dot and tunnel barrier using 
Eq. (fOl) 

Qr = —{Qd + Qt) ■ (14) 

We recall that our calculation applies only for the linear 
response regime and thus we keep only the leading order 
contributions in the fields V (t) and U (t). We employ the 
nonequilibrium Keldysh Green’s function formalism for 
the calculation of the energy change rates in timei^ We 
start with the time derivative of Rp first, 

J^{nD{t))= {dtemnAt)]{dlAt)dnAt)) 

m,n,(T 

+ ^ 'dmnait)dt{d\^„{t)dna{t)) . (15) 

m,n,cr 


dt 


( 12 ) 
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The first term on the right hand side denotes the power 
developed by the AC source and does not contribute to 
the heat flow. The second term represents the energy 
flux going into the dot which can be written in terms of 
the dot Green’s function 


Qoit) = -fTr . (16) 

Hereafter, bold face symbols denote matrices such that 
their (m, n) components are 


[^da{t)] 
[Gda,da(j'y ^0] 


m,n 

Tn,n 


^mn<7 


(17) 


The trace means summations over energy levels and spin 
indices, i.e., Tr = n'^a- Secondly, from the defini¬ 

tion of it can be shown that 


{-Ht) = -iTi t)t, + t:G< t)] , (18) 

and thus the time variation of energy stored or relaxed 
at tunneling barrier Qrit) = dtiT-Lr) is given by 

Qrit) = -idtTr [t:G< ,,,(t, t) + G<,,,,(t, t)t,] . (19) 

Here, the reservoir-dot and dot-reservoir Green’s func¬ 
tions are given, respectively, by 


{Gda,cait^t')]n,k = i{cl„{t')dn^ (t)) , 


and and Tr = Y.k,n,a- 

Using the standard technique of equation of motion, 
the reservoir-dot Green’s function can now 

be recast in terms of solely the dot Green’s functions as 
follows. 


being f(ti — ^ 2 ) the Fermi-Dirac distribution function in 
the time domain. Similarly, the dot-reservoir Green’s 
function is given by 

G,V(t,t)G = J ^ [GL,,.(i,U)S<(ti,t) 

+ Gt,dAt,h)Kitut)]. (25) 

Employing Eqs. (EH and EH, we find Qt is also ex¬ 
pressed only in terms of the dot Green’s functions. 

Once the heat energy rates for the dot and the tun¬ 
neling parts (Qd, and Qt) are expressed in terms of the 
dot Green’s function, the next step consists of comput¬ 
ing such Green’s functions in the presence of the time- 
dependent AG signal. For the dot’s retarded/advanced 
Green’s function, we have^i 

GtdaM = (26) 

where = J^,{dti/h) eU{ti) and - ^ 2 ) 

denotes the dot’s retarded/advanced Green’s function in 
equilibrium, i.e., in the absence of AG signal. Finally, the 
dot’s lesser Green’s function is obtained by means of 

Gt,dAt,t') = 

GS,,,,(t,U)S<(ti,t2)Gl,,,(t2,t')- 

(27) 

Inserting Eqs. (I21l) - (l27p into Eqs. (1161) and E3 , it can 
be shown, after some cumbersome algebra, that to linear 
order in V and U the energy change rates in frequency 
domain become 


t:G 


< 

C<T,d <7 


(t,t) = J ^ [S;(t,ti)G,V(U,t) 

+ S<(t,H)G^,,,,(U,t)]. (21) 


Qoiuj) = zewTr| J ^ edcrA{e, ftw) [U(fcu) - t/(feu)] | 


(28a) 


Here, we have defined 

^:ih,h) = -'-hT6{h-t2), ^ 22 ) 

S<(ti,t2) = ire-*'^'^(‘i’‘^V(U-t2) 

as the retarded and lesser self-energiesii that contain the 
time-dependent fields. We have assumed a momentum- 
independent tunneling amplitude tnka- = tna- leading to 
hybridization strength given by Pot^^^tno- 

with po = 1/211 the density of states of the reservoir 
and D the reservoir bandwidth. In such self-energies, it 
is defined 

/•‘i dt 

(I)v{ti,t 2 ) = J —eV{t), (23) 

as the time-dependent phase due to the AG external po¬ 
tential and 


Qt{uj) = ieu;Tr| 


— (Sw -b 2e - 2ed^)A{€, fiw) 

x[V{huj)-U{nuj)]'^ , (28b) 


where V{huj), and U{huj) are the Fourier transforms of 
U(t), and U{t) respectively. Here, 


A{e, huj) = T{huj, e)J-{e, huj) 


with 


T{haj, e) 
J-'(e, hbj) 


Gr.L(e + ft^)rG“/^L(e) 

f{t + fiw) - /(e) 


(29) 


(30) 


The dot’s retarded/advanced Green’s function in equilib¬ 
rium and in frequency domain is given by 


f{ti-h)= I , (24) 


f^r/a 

^d<7,da 


(e) 


1 


e — edcr ± fr/2 


(31) 
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with 

[^(icr]m,n — ^m,n I ^mcr “t“ ^ ^ I 

y l,s{^m,<7) j 

- 2i?c(dLd„<,). (32) 

Although Eq. (IHTl) has the same structure with its nonin¬ 
teracting retarded/advanced Green’s function, the charg¬ 
ing energy is included in the denominator such that it 
can describe the Coulomb blockade effect as mentioned 
above. To completely determine the dot Green’s func¬ 
tion, the diagonal and off-diagonal dot occupations need 
to be self-consistently calculated using 

i^hs^ms) =J ■ (33) 


immediately generalized to our continuum model by just 
keeping the factor one-half in front of the tunneling en¬ 
ergy flow. 

After these considerations, the formulations for the the 
dot and reservoir frequency dependent heat currents in 
the linear response regime read 

Jd = - Jr = *ewTr| j ^ -|- A{e, ftu) 

X \y{hu) — C/(fiw)]|. (36) 

Importantly, this formula is the central finding of our 
work. 


D. Electrothermal admittance 


Finally, from Eqs. (l28ll and o, we obtain the expres¬ 
sion for the energy change rate at the reservoir 

Qit(w) = -{Qd + Qt) = -jea;Tr| J ^ 

{huj + 2e — eda) A{e, huj) [V(hw) — U{hu)] | . (34) 

Importantly, neither Q_r(w) nor Q_d(u;) have a well de¬ 
fined parity when the AC frequency is reversed, and 
therefore, within linear response theory, these two magni¬ 
tudes do not represent physical quantities [see Sec. im. 
Based on these observations, the reservoir and dot fre¬ 
quency dependent heat current expressions must be thus 
modified to exhibit a proper parity property when the 
AC frequency is reversed. We find the expressions 

Jr = Qr{^) + 2 Jd = Qd(w) -I- (35) 


In the presence of a time-dependent driving force, it 
is quite general to characterize the transport using the 
concept of admittance. The complex electrothermal ad¬ 
mittance is defined as 

= W)- 

Notice that Eq. (1571) which can be obtained from Eq. (1551) 
contains the unknown function Uifiw). Therefore, for a 
complete characterization of the complex electrothermal 
admittance, we need hrst to determine the internal po¬ 
tential U{fujj). For such purpose, we note that the dis¬ 
placement current Xu can be featured by a capacitance C 
in a simple model [here, we consider the situation shown 
in Fig. [fl (a)] 

Xd{uj) = -iujCU{nuj). (38) 


satisfy the parity property. Remarkably, the choice of the 
factor ^Qt(w) is unique in order to ensure a well defined 
parity property in both Jr, and Jd. 

Formally, these expressions agree with their non¬ 
interacting counterpart for the time dependent heat 
currentsi^ However, it should be noted that our theo¬ 
retical analysis (i) goe s beyond low AC frequencies in 
contrast with Ref. in which the heat rate was ob¬ 
tained up to second order in the AC frequency using the 
Floquet theory, (ii) includes the effect of Coulomb block¬ 
ade, and (Hi) is applicable to multi-orbital conductors in 
contrast to previous time heat formulations. 

Remarkably, our AC heat formula contains photon- 
assisted tunneling events only possible for sufficiently 
high AC frequencies. Note that a similar definition for 
the time-dependent heat, but applicable to a spin chain 
model (using a tight-binding model), was proposed in 
Ref. [^. In such work, the heat flux connecting two 
sites, i and z ± 1, was incorporated to the general heat 
flow expression with a 1/2 factor in close analogy to 
Eq. (1351) . However, caution is needed when this com¬ 
parison is made. The results for a chain of sites are not 


Due to current conservation, the displacement current 
is equal to the tunneling current for a quantum 

capacitor 


= g{uj)[V{hw) - U{huj)], (39) 

where [see Ref. [13 for its explicit derivation] 

g(a;) = ze^wTrjy ^A(e,/ia;)| . (40) 


The internal potential is obtained when we impose cur¬ 
rent conservation Xu =It, then U{huj) reads 


u{huj) 


g{uj)V{fu^) 
—iivC + g{oj) 


(41) 


Inserting Eq. (1511) into Eq. (1551) completely characterizes 
the linear response of the heat current to a time depen¬ 
dent voltage, i.e., the electrothermal admittance 


= m{uj) 


—iuiC 


[-iujC + g{uj)] ’ 


(42) 





6 


with 


m{uj) 

Similarly, 


= iewTr Iy ^ | . (43) 

the thermoelectrical admittance is defined as 

C(u^) = (44) 

T{uj) 


with T(uj) being the Fourier transform of a time modu¬ 
lated temperature T{t). Remarkably, the thermoelectri¬ 
cal and the electrothermal admittances are reciprocally 
related due to the microreversibility principle 


M{uj) = TC{uj). (45) 


Notice that microreversibility principle only holds at lin¬ 
ear order in V{t), and U{t). Finally, a second order ex¬ 
pansion in the AC frequency for m{uj) 

m{ui) = —iojCm + , (46) 


allows us to obtain the RC electrothermal parameters. 
For comparison, we also calculate the electrical admit¬ 
tance defined by 


Giu) 


Xj- — iu)C 

Vihjj) ^ + 


(47) 


The corresponding second order expansion of g{uj) in fre¬ 
quency 


g{uj) = -iujCg + uj^C^Rg, (48) 

always yields positive RC parameters, i.e., It is always 
delayed with respect to R(w). This is in clear contrast 
to the electrothermal admittance case in which both Cm 
and Rm can be either positive or negative. The heat 
flow response is elapsed or delayed with respect to the 
electrical signal depending on the system parameters. 
Note however that the product of both quantities, the 
RC time, is kept always positive as expected. Similar 
results were obtained in Ref. [2^. 


III. RESULTS 

Previously, we have derived the formal expressions for 
the reservoir heat flow and the corresponding electrother¬ 
mal admittance. Now, we apply these formulas to the 
case of an interacting multi-orbital conductor. However, 
before addressing the multi-orbital case, for illustration 
purpose we first investigate a single orbital interacting 
quantum capacitor under the influence of an AC poten¬ 
tial. Later on, we will analyze a prototype of multi¬ 
orbital conductor, a carbon nanotube quantum dot at¬ 
tached to a single reservoir. In both cases, we show re¬ 
sults for the electrical and electrothermal transport. 



Ej/r 


FIG. 3: Charge relaxation resistance for a single orbital qnan- 
tum capacitor versus the Fermi energy. Parameters: ed ~ 0, 
Az = O.SF, 2Ec = lOF, and ksT = 0.04F. 



FIG. 4: (a) Real and (b) imaginary parts of gigj) for a single 
orbital quantum capacitor as a function of the AC frequency. 
Photon-assisted excitations occur at the transition rate ener¬ 
gies ±huj = \Ef — Edo-1 when ksT <C P. Parameters: td = 0, 
Az = 0.625r, 2Ec = lOP, Ep = 0, and ksT = 0.04F. 


A. Single orbital quantum capacitor 

We consider an interacting quantum dot with just one 
orbital contacted to a single reservoir. The reservoir has 
an electrical potential that oscillates in time. In accor¬ 
dance with our previous theoretical considerations, the 
prefactor g{ijj) of the electrical admittance I/(a;) can be 
expanded in powers of the AC frequency up to second 
order as shown in Eq. (1481) . A very well known result 
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FIG. 5: (a) Real and (b) imaginary parts of m(cj) for a single 
orbital quantum capacitor versus the AC frequency. Photon- 
assisted excitations occur at the resonant conditions. Param¬ 
eters: td = 0, Az = 0.625r, 2Ec = lOP, Ep = 0, and 
O.OdP. 


establishes that Rg = hj{2Qe^) becomes universal be¬ 
ing Q the number of transport channels. Our example 
as shown in Fig. |3] considers the presence of a Zeeman 
field to explicitly break spin degeneracy. Rg takes the 
value of /i/2e^ when just one of the two spin-resolved 
levels (ed-f/j. = Ep ± Azl2) participates in transport, 
whereas it becomes h/4e^ when the opposite spin channel 
also contributes. Clearly, the two main peaks observed 
in Rg are separated roughly by 2Ec the charging energy 
{2Ec 2> Az). Our results for Rg indicate the presence of 
charging effects and the Coulomb blockade phenomenon. 

More importantly, our major interest resides in the be¬ 
havior of the electrical and electrothermal admittances at 
arbitrary AC frequencies. Figured shows the real [Fig.0] 
(a)] and imaginary [Fig. 0] (b)] parts of the prefactor 
g{uj) in the expression for the electrical admittance f/(w) 
[see Eq. gzl)]. We observe that ^eg{oj) = 5Reg(—w) and 
'^mg(uj) = —3mg(—w). The factor—fwC'/(—*wC'-|-g(w)) 
in Eq. dUl) then does not change the parity property of 
Q{uj) with respect to uj such that for simplicity we con¬ 
sider g{uj). Excitations occur when the AC frequency 
matches with the resonant condition zthuj ss \eda — Ep\. 
We recall that Cda is the spin-dependent dot energy level 
renormalized by electron-electron interactions according 
to 

€da = Cd + Az + 2Ec{d^ddd) ■ (49) 

Eor the parameters used in Eig. 01 we obtain 'eda ~ 
O.SF, 2.IF which agrees with the resonant behavior found 
in g{uj) with peaks at Ahcu « O.SF, and Ahu ^ 2.IF. 

Similar features are observed in the electrothermal ad¬ 
mittance shown in Fig. [5] Interestingly, the imaginary 
part of m(uj) for w > 0 takes either positive or negative 
values by tunning the AC frequency, which indicates that 
time heat current can be either delayed or elapsed with 


respect to the AC electrical time-dependent signal. 

Now, we discuss the parity property of the response 
functions g(aj) and m(uj) [thus 0(uj) and A4(a;)]. We 
write g(uj) in the form 

g(uj) = Sfteg(Lj) + iQmg(uj) , (50) 

and express the real/imaginary part as 
1 1 /■°“ 

^eg(u;) = - Ig(uj) + g*(uj)] = -J dt [g{t) + g{-t)] , 

(51) 

^mg{u) = — [ 5 r(w) - g*{uj)] = J [^(i) - 9{-t)] ■ 

(52) 

Here, we used the fact that the response function g{t) 
must be real to have a real expectation value for the 
current Ipit)- From Eq. (l52l) . it is quite easy to show 
that 

1 r°° 

^mg{-u}) = Yi J * - git)] , ( 53 ) 

which implies ^mgiui) = —5mg(—w). Using a simi¬ 
lar line of reasoning, we can also prove that ^eg{u;) = 
5fteg(—w). This argument also works for m(a;). It is wor¬ 
thy to notice that this parity property comes from the 
fact that we have included the contribution due to the 
tunnel Hamiltonian in our definition for the heat flow 
as shown in Eq. (I35|) . Eurthermore, this parity argu¬ 
ment goes beyond the simple site partitioning scheme 
explained in Ref. [s^. 

B. Multi-orbital quantum capacitor 

We now investigate a single reservoir carbon nanotube 
quantum dot as an example of a multi-orbital interact¬ 
ing conductor. We regard the nanotube quantum dot as 
a localized single particle level described by two quan¬ 
tum numbers, the orbital quantum number r associated 
to clockwise (r = -|-1) and anti-clockwise (r = —1) semi- 
classical orbits along the nanowire circumference (related 
with the AT-valley degeneracy in graphene) and the spin 
degree of freedom a. In the presence of magnetic field 
along the nanotube axis, the dot energy level splits in 
the spin sector by the Zeeman field Az and in the or¬ 
bital sector by an amount Aorb that depends on the nan¬ 
otube radius^i Besides, due to the nanowire curvature a 
non-negligible spin-orbit interaction is present, yielding a 
Kramers splitting of magnitude Ago. All together yields 
the following carbon nanotube dot energy level: 

Cdar = Cd + O-A^ -I-rAorb + CTtAso. (54) 

In the following, we show results that correspond to the 
realistic parameters!^ = 0, A^ = 0.625r, Aorb = 
5 A 2 (orbital magnetic moments can be 5 — 20 times 











larger than its spin counterpart^), and Ago = 0.5Az- 
The charging energy is 2Ec = 2r which lies in the strong 
interacting regime. We restrict ourselves to the low tem¬ 
perature regime ksT = 0.05r. 

Ep/T Ep/T 



FIG. 6: (a) Electrical capacitance Cg and (c) electrical re¬ 
laxation resistance Rg versus the Fermi energy Ep. (b) Elec¬ 
trothermal capacitance Cm and (d) electrothermal relaxation 
resistance Rm versus the Fermi energy Ep- Nanotube pa¬ 
rameters: ed = 0, Az = 0.625r, Aorb = 5Az, Aso = 0.5Az, 
2Ec = 2r, and kpR = O.OSF. 


two consecutive resonance points. The sign inversion in 
Rm happens when the derivative of the carbon nanotube 
quantum dot density of states vanishes.— Remarkably, 
both Cm and Rm diverge around the resonance points, 
behavior that is washed out by enhancing temperature. 
(Similar results were obtained for the weak interacting 
limit, when Ec -C T see Ref. [2^ for details). 



FIG. 7: (a) Real and (b) imaginary parts of g{oj) versus the 

AC frequency uo. Parameters: ea = 0, Az = 0.625r, Aorb = 
5Az, Ago = 0.5Az, 2Ec = 2r, Ep = 0, and ksT = O.OSF. 


First, we show our results for the electrical and elec¬ 
trothermal RC parameters in Fig. [51 In the low tem¬ 
perature limit {ksT <g; F) the RC parameters exhibit 
universal values.— We recall that in the low frequency 
regime the RC parameters characterize the electrical and 
electrothermal admittances [see Eqs. (gni) and (gSl)]. Cg 
and Rg are displayed in Fig. [5] (a) and (c) when Ep is 
varied. The electrical capacitance Cg shows oscillations 
which peaks at the positions located roughly at 

€d- Az - Aorb + Ago « -3.4r, (55) 

ed + Az — Aorb — Ago -I- 2Ec « —O.SF, 

ed — Az + Aorb — Ago -I- 2{2Ec) ~ 6.2r, 

ed + Az + Aorb + Ago -l- 3(2Ec) ~ lOF. 

As expected, at each nanotube level position, Rg takes 
the value of /i/2e^, while in the middle of two consecu¬ 
tive resonances (when two resonances contribute to Rg) 
it diminishes to half of this value. Whereas Rg and Cg 
display always positive values, the electrothermal capaci¬ 
tance Cm and resistance Rm can become positive or neg¬ 
ative when Ep is tuned. This is shown in Fig. |n](b) and 
(d). Indeed, Cm changes sign whenever the Fermi en¬ 
ergy matches with any of the nanotube resonances. Heat 
current becomes delayed or elapsed with respect to the 
AC signal depending on the Fermi energy position. The 
electrothermal resistance Rm modihes its sign not only 
at the points when Ep coincides with the nanotube res¬ 
onances [see Eq. (1551) ] but also when the electron-hole 
symmetry point occurs, just at the midpoint between 


We now discuss the case of arbitrary AC frequencies 
and analyze the prefactor g{u}) of the electrical admit¬ 
tance in Fig. 0 The real and imaginary parts of g{ui) 
versus the AC frequency are depicted in Fig. [7| (a) and 
(b). As in the single orbital quantum capacitor, we 
observe that ^mg{u}) has odd parity with respect to 
Lo, while 'Steg{Lu) is an even function of w as a conse¬ 
quence of being response functions of a real perturbat- 
ing force. ^mg{uj) accounts for the dissipative part of 
the electrical conduction with resonances roughly located 
at zLhuj Ki \Ep — Qo-tI- For the HF approximation, 
these resonances coincide with the dot level positions that 
are renormalized by interactions according to Eq. dloj. 
These renormalized level positions (when Ep = 0) are at 
edar ~ — 1.4r, —O.SF, 6.IF, 8.OF leading to the observed 
resonances in g{u;). The resonant behavior of g{uj) re¬ 
flects the photon-assisted tunneling processes in which 
transport through the nanotube occurs by absorpting or 
emitting single photons. Furthermore, these resonances 
are broaden mainly by F at very low temperatures. The 
3?eg(w) corresponds to the reactive part of the electrical 
conduction and has also a similar resonant structure. 

We turn into the analysis of the electrothermal admit¬ 
tance Mioj) which is given in Eq. (H^ . As before, it is 
convenient to examine m(uj) defined in Eq. (1431) and we 
thus plot it as a function of the AC frequency in Fig. |S| 
First, we observe m{uj) displays the same parity property 
as g{ijj) which can be explained in the same way. Fur¬ 
thermore, the real and imaginary parts of m^uj) are also 
characterized by features located at ±huj ^ \Ep — edar]- 
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FIG. 8: (a) Real and (b) imaginary parts of m{ij) versus the 
AC frequency lj. Parameters: Sd ~ 0, Az = 0.625r, Aorb = 
5Az, A,o = 0.5Az, 2Ec = 2r, Ef = 0, and ksT = O.OSF. 


Remarkably, the observed peaks indicate single photon 
absorption and emission processes for the transport of 
heat. It is worthy to emphasize that our approach is 
able to capture such photon-assisted processes in con¬ 
trast with other calculations restricted to AC frequencies 
smaller or similar to the tunnel coupling F: the other 
calculations assume the AC frequency is small such that 
only a second order expansion of m(aj) in powers of oj 
(characterized by the pairs Cm, Rm) can be justified. 
However, since we do not have such a restriction for the 
applied AC frequency, our results clearly exhibit photon- 
assisted processes for the heat transport. 

Remarkably, m(uj) (either its real or imaginary part) 
takes positive or negative values depending on the AC 
frequency regime. We observe 3fiem(w) < 0 for low and 
moderate w, whereas for large AC frequencies its sign is 


reversed. We stress the importance of such result for the 
functionality of quantum circuits, in which one could ma¬ 
nipulate the sign of the heat flow spectrum by properly 
tunning the AC frequency. 


IV. CONCLUSIONS 

In closing, we have investigated the heat current spec¬ 
trum in the linear response regime for an interacting con¬ 
ductor coupled to a single reservoir and modulated by 
an electrical AC signal. Our results, valid for arbitrary 
AC frequencies, show that the heat current expression 
for the reservoir needs to consider the heat stored or re¬ 
laxed at the barrier. We illustrate our findings with two 
prototype of interacting conductors, namely, a single or¬ 
bital quantum dot and a multi-orbital conductor~a car¬ 
bon nanotube quantum dot coupled to a single reservoir. 
We deal with the strong interacting limit where Coulomb 
blockade phenomena applies. We highlight that (i) elec¬ 
trical and heat transport displays photon-assisted trans¬ 
port features, and (ii) the electrothermal admittance can 
be positive or negative and the sign can be chosen by ad¬ 
justing properly the AC frequency. This is an important 
issue for engineering nanoelectronic circuits with optimal 
heat dissipation performances. 
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